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Abstract 

"^rhe  feasibility  of  a  red  singlet  oxygen  laser  based  on  production  of  (b) 
state  singlet  oxygen  by  the  electronic  energy  pooling  reaction  of  (a)  state 
singlet  oxygen  was  Investigated.  A  microwave  discharge- flow  system  for 
studying  the  quenching  kinetics  of  (b)  state  obygen  was  constructed  with  a 
large  diameter  viewing  region  so  as  to  ensure  that  heterogeneous  quenching  of 
both  (a)  and  (b)  state  oxygen  was  unimportant. 

An  upper  bound  for  quenching  of  (b)  state  oxygen  by  (a)  state  oxygen  was 
measured  and  was  not  found  to  be  extremely  large  compared  to  the  rate  constant 
for  quenching  of  (b)  state  oxygen  by  ground  state  oxygen. 

The  rate  constants  for  quenching  of  both  (b)  state  oxygen  and  (a)  state 
oxygen  by  ground  state  oxygen  was  measured.  Both  rate  constants  agreed  well 
with  those  measured  by  previous  workers. 

The  results  suggest  that  the  extent  of  the  effective  population  Inversion 
of  (b)  state  oxygen  produced  energy  pooling  will  not  be  suppressed  by  the  high 
densities  of  (a)  state  oxygen  Initially  present.  The  measurements  reported 
here  support  the  feasibility  of  the  red  singlet  oxygen  (Red  SOX)  laser  and 


further  research  and  development  toward  this  goal  Is  recommended. 
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A.  Introduction 

Chemical  lasers  are  among  the  devices  being  considered  for  delivery  of 
high-power  directed  energy.  The  high  energy  storage  density  potentially 
available  with  chemicals  makes  chemical  lasers  particularly  attractive  for 
airborne  missions,  where  weight  considerations  are  Important.  However, 
chemical  systems  presently  under  evaluation  suffer  from  a  number  of  drawbacks 
which  are  expected  to  be  greatly  reduced  In  the  Singlet  Oxygen  (SOX)  laser 
which  Is  the  subject  of  this  proposal. 

Chemical  lasers  fall  Into  two  principal  groups,  depending  on  whether  the 
spectroscopic  transition  Involves  vibrational  excitation  of  a  molecule  or 
electronic  excitation  of  an  atom  or  molecule.  Examples  of  the  former  laser 
group  are  the  HF  and  DF  lasers,  pumped  by  reacting  fluorine  or  a  fluorine- 
containing  compound  with  H2  or  02.  These  lasers  are  plagued  with  problems  of 
fluid  dynamics  associated  with  high  energy  release,  extremely  rapid  rates  for 
reaction,  and  comparatively  fast  colllslonal  deactivation  rates.  In  addition, 
the  long  wavelengths  (2.4  -  3.8  urn)  present  problems  of  atmospheric  trans- 
mission  and  low  brightness,  the  latter  of  which  scales  with  x  . 

The  only  known  example  of  a  scalable  electronic  transition  laser  Is  the 
chemical  oxygen-iodine  laser  (COIL) . In  which  electronic  excitation  of 
metastable  02(1a)  molecules,  produced  by  chemical  reaction  of  chlorine  with 
basic  hydrogen  peroxide.  Is  transferred  to  atomic  Iodine.  The  laser  operated 
on  the  Iodine  transition  at  1.3  urn.  While  the  wavelength  of  the  COIL  laser 
permits  over  eight  times  the  brightness  of  a  DF  laser  operating  at  the  same 
power  level,  frequency  doubling  and  tripling  are  still  being  considered  to 
further  Improve  brightness.  Other  Important  issues  In  COIL  are:  the  need  for 
supersonic  flows  to  Improve  medium  homogeneity  and  counteract  high  quenching 
rates,  the  uncertainty  about  the  mechanism  of  I2  dissociation  and  thus  the 


1 


scalability  to  higher  densities,  problems  In  obtaining  uniform  mixing  on  a 
short  time-scale,  and  the  need  for  generators  of  02(*a)  that  will  operate  at 
low  gravity  and  variable  acceleration. 

There  Is  thus  an  Impetus  to  look  for  new  chemical  lasers  that  might  have 
shorter  wavelengths  and  fewer  fluid  dynamic  and  kinetic  rate  limitations  than 
current  systems  under  consideration.  A  laser  based  on  conversion  of 
chemically-produced  02(1a)  to  02(1e)  Is  proposed  as  a  viable  candidate.  The 
lasing  wavelength  of  0.76  urn  could  provide  three  times  the  brightness  of  a 
COIL  laser.  The  expected  slowness  of  quenching  reactions  as  well  as  primary 
production  would  greatly  reduce  hydrodynamic  and  kinetic  mixing  constraints 
and  provide  for  large,  uniform  gain  regions. 

In  the  past,4*5  attempts  have  been  made  to  obtain  direct  lasing 
of  02(1a),  motivated  by  a  desire  to  circumvent  complications  that  arise  from 
mixing  of  Iodine  In  a  COIL  laser.  These  attempts  have  been  unsuccessful 
because  of  the  low  cross-section  for  stimulated  emission  and  thus,  low  gain  of 
an  02(*a)  laser.  This  difficulty  could  be  greatly  reduced  In  a  laser 
utilizing  02(1l) 


O^e)  *  02(3e)  +  hv 

because  the  cross-section,  and  thus  the  small  signal  gain  coefficient  for  the 
Doppler-broadened  transition.  Is  over  50  times  greater  than  for  the  transition 
In  02(1a). 

The  source  of  02(*e)  for  a  laser  would  probably  be  the  energy  pooling 
reaction 


(1) 


2 


V.N  Va 


As  set  forth  In  the  Phase  I  Statement  of  Work,  the  primary  goal  of  this 
Investigation  has  been  to  measure  the  rate  constant  for  quenching 
of  02(1e)  by  02  (1a); 

02  (lt)  *  02  (M  ---*  02  (h)  *  02(1a)  (2a) 

k 

™  02  (h)  +  02(3e)  (2b) 

Reactions  (2a)  and  (2b)  are  potentially  Important  sources  of  deactivation 
of  02(1r)  under  the  Initial  conditions  of  high  02(*a)  density.  The  production 
of  an  electronic  population  Inversion  In  02  requires  that  (2a)  and  (2b)  be 
slow  reactions  relative  to  Reaction  (1).  The  aim  of  Phase  I  of  this  project 
has  been  to  quantitatively  assess  the  Importance  of  Reactions  (2a)  and  (2b)  In 
gaseous  oxygen  systems  by  obtaining  an  estimate  to  k£2. 

We  next  discuss  the  results  of  our  Phase  I  Investigation  and  their 
significance  to  the  Phase  II  effort.  In  Section  B  we  will  discuss  the 
experimental  apparatus  and  procedures  used  to  perform  the  measurements.  In 
Section  C  we  give  a  kinetic  analysis  for  the  steady  state  and  time  dependent 
density  of  02(*a)  and  02(*e).  The  experimental  results  are  presented  and 
discussed  in  Section  D. 

B.  Experimental  Apparatus  and  Procedures 

The  measurements  have  been  performed  using  a  discharge  flow  apparatus 


with  a  large  (60  cm)  diameter  viewing  region,  a  schematic  of  which  Is  shown  In 
Figure  I.  Flow  rates  were  measured  using  a  calibrated  mass  flow  meter  (Scott) 
and  pressure  was  measured  to  within  0.02  Torr  by  a  capacitance  manometer 


Singlet  oxygen  was  produced  by  a  2.45  GHz  microwave  discharge  through 
four  parallel  3.9  an  l.d.  quartz  tubes.  Surrounding  each  quartz  tube  was  a 
microwave  cavity  containing  three  magnetrons  (Toshiba  2M205  and  2M207)  driven 
by  a  208  V  three  phase  circuit  such  that  each  magnetron  within  a  group  was 
driven  120°  out  of  phase  with  respect  to  the  other  two.  By  monitoring  the 
02  (*a)  and  02  (M  density  downstream  of  the  discharge  we  were  able  to 
verify  that  this  discharge  arrangement  produced  steady  densities  of  02(1a) 
and  02(1r)  with  no  evidence  of  modulation  at  60  Hz. 

The  density  of  02(*a)  was  measured  with  a  germanium  photodiode  (Judson 
J-16D)  cooled  to  77  K  and  calibrated  against  a  standard  tungsten  lamp.  An 
Interference  filter  with  approximately  90  nm  bandpass  Isolated  the 
1.27  v  (a  ♦  t)  fundamental  emission  of  02(1a).  The  current  output  of  the 
photodiode  was  measured  using  a  Kelthley  410A  plcoammeter. 

Both  detectors  were  mounted  on  carriers  which  moved  along  an  optical  rail 
parallel  to  the  vessel  axis.  The  02  (lr)  density  was  measured  with  a 
Hamamatsu  R1104  photomultiplier  tube  calibrated  against  the  same  standard  lamp 
and  equipped  with  a  5.5  nm  bandpass  Interference  filter  centered  at  760  nm 
which  Isolated  the  b+X  transition.  The  PMT  signal  was  observed  on  an 
oscilloscope  with  1.00  Mq  Input  Impedance.  Each  detector  was  equipped  with 
Imaging  optics  so  that  the  center  of  the  viewing  vessel  and  the  detector  faces 
were  conjugate  points.  This  arrangement  provided  for  uniform  weighting  of 
densities  at  each  radial  position  In  the  vessel,  l.e.  the  detected  signal 
for  species  1  Is  given  by 

R 

Si  -  C,  f  0,(r)dr  (3) 

•R 

In  which  p^(r)  is  the  density  of  species  1  and  Is  a  calibration  factor 


relating  density  of  species  1  to  signal  strength.  R  Is  the  tube  radius.  Note 
however,  that  If  p^(r)  constant,  then  Equation  (3)  does  not  properly  weight 
the  density  p^(r)  to  give  a  signal  proportional  to  the  average  density  over 
the  tube  cross-section.  The  latter  would  require  an  additional  factor  of  r 
under  the  Integral  In  order  to  weight  regions  near  the  vessel  wall 


Si  -  ci  J  rPj(r)dr  .  (4) 

0 

2  2  2 

For  parabolic  distributions  of  the  type  p(r)  *  pQ(r  +  const)  or  p(r)  =  p  (R  -r  ), 
which  are  expected  when  species  are  produced  or  destroyed  at  the  wall,  the 
difference  between  Equations  (3)  and  Eq.  (4)  can  be  large.  Measurements  were 
made  to  determine  whether  densities  were  strongly  dependent  on  r,  and  these 
measurements  will  be  reported  In  Section  30. 


For  the  present  Investigation  it  was  important  to  verify  that  quenching 
of  O^a)  and  O^1!)  by  species  other  than  0£  was  not  Important.  The  species 
which  posed  the  greatest  potential  for  complicating  our  quenching  kinetics 
results  were  NO2,  O3,  0,  and  H2O.  Blmolecular  rate  constants  for  quenching 
and  other  reactions  involving  O^a)  and  Og^x)  are  given  in  Table  1.  We  have 
evaluated  the  effects  of  these  species  as  follows. 


HI) 

_ 2.  Nitrogen  dioxide  Is  formed  In  oxygen  discharges  when  N2  is  present  as  an 

Impurity  In  commercial  oxygen.  NO2  Is  a  known  quencher  of  02(1a)  and  02(*x), 
with  quenching  rate  constants  of  5  x  10“ 18  and  2  x  10“14  cm3  molec-1  sec-1, 
respectively  (See  Table  1).  Assuming  the  extreme  case  that  all  of  the  N2 
present  In  a  1  To rr  flow  of  oxygen  (<  10  p.p.m.)  is  converted  to  NO2,  the 
pseudo-first  order  rate  constants  for  quenching  02(1a)  and  02(1r)  are  6.4  x  10“6 


and  3.2  x  10"2  sec"*,  respectively.  The  corresponding  pseudo-first  order  rate 
constants  for  quenching  of  02(*a)  and  02(*z)  by  ground  state  oxygen  at  1  Torr 
are  5  x  10"*  and  1.4  sec"*,  respectively.  We  conclude  that  quenching  by  N02 
Is  not  significant  In  this  experiment. 


S***.  In  low  pressure  oxygen  discharges,  oxygen  atoms,  molecular  oxygen  and 
ozone  participate  In  a  coupled  kinetic  scheme.  The  rate  constants  for 
quenching  of  02(*a)  and  02(*z)  by  O3  and  0  are  given  In  Table  1.  The  dominant 
homogeneous  reactions  occurring  downstream  of  the  discharge  governing  the 
production  and  destruction  of  these  species  can  be  summarized  by  the  following 
kinetic  scheme. 


0  +  o2  +  o2  o3  +  o2 


1  kg 

03  +  02  (Aa)  -S  0  +  02  +  02 


o3  +  0  — »  Op  +  02 


The  rate  equations  describing  the  disappearance  of  0  and  O3  are 
“dlOJ  2  1 

— ft*-  *  -kg  [01  [02r  +  kg  [03I {02aa|  +  k7[03l [0|  (8) 

=$*-  *  k5(0][02J2  -  kg  103!(02*a]  +  k7  [031(01  (9) 

Since  02  (*z)  and  02  (*a)  are  the  dominant  species  present  In  the  discharge. 
Equations  (8)  and  (9)  may  be  rewritten  In  pseudo-first  order  form 


-d(03J 


-  k5  [0]  +  kg  [03l  +  k7  [03] [01 


-g^=k5  [01  -  k6  I031  ♦  k7  IO3IIO]  (11) 

■*  9  1 

where  kg  *  kg  [ 0^ 1  and  kg  =  kg  (02  a!.  Using  the  literature  values  (Table  1) 
for  kg  and  kg,  and  using  typical  values  for  O2  and  O^1*)  of  1  Torr  and  0.1 
Torr  respectively,  we  find  that  kg  =  0.4  see"*  and  kg  =  11  sec"1. 

Reactions  (5)  and  (6)  do  not  result  in  net  removal  of  O3  and  0.  In  the 

gas  phase  removal  is  governed  by  Reaction  (7),  that  is,  the  direct  reaction  of 
O3  and  0,  with  k7  =  6.2  x  10"15  cm"1  molec"1  sec-1.  Assuming  equal  initial 
densities  of  O3  and  0  (say  1014  cm-^)  the  decay  will  be  second  order  with  the 

first  half  life  (k7(0])_1  =  1.6  sec.  The  decay  times  of  02(1a)  and  02(1r) 

are  comparable  to  or  shorter  than  this  at  O3  and  concentration  of  1014  cm“^. 

We  conclude  that  gas  phase  reactions  will  not  quantitatively  remove  O3  or 
oxygen  atoms  sufficiently  rapidly  to  prevent  their  quenching  of  02(1a) 
and  02(1r). 

In  addition  to  removal  of  O3  and  0  by  homogeneous  reaction,  oxygen  atoms 
recombine  on  the  walls  of  the  discharge  tube.^2  Formation  of  a  post-discharge 
HgO  film  on  the  inner  surfaces  of  the  quartz  tubes  by  addition  of  Hg  to  the  02 
flow  prior  to  discharge  Is  known7,24,2®  to  increase  the  rate  of  0  atom 
recombination.  During  the  course  of  this  study  the  effect  of  operating  the 
discharge  with  and  without  the  HgO  film  was  Investigated.  It  was  found  that 
with  the  HgO  film  In  place,  the  steady  state  density  of  02(1a)  was  not 
significantly  different  from  that  measured  with  no  film,  however  the  density 
of  02(1i)  was  an  order  of  magnitude  larger.  Given  the  large  rate  constants 
for  quenching  of  02(1t)  by  03  and  0  we  conclude  that  the  HgO  film  probably 
significantly  reduces  O3  and  0  density,  but  the  evidence  Is  Indirect,  and  we 
did  not  have  the  capability  to  establish  the  presence  of  these  species  nor  to 
determine  whether  removal  of  these  species  Is  nearly  complete  in  the  presence 
of  the  HgO  film. 


-12  ~_3 
cm 


HO  1 

2 .  The  rate  constant  for  quenching  of  02(  z)  by  H20  is  6  x  10 

molec"1  sec"1  (see  Table  1).  For  this  reason  we  have  taken  care  to  exclude 

traces  of  water  from  the  flow  system  by  passing  02  through  a  stainless  steel 

coil,  20 0  cm  In  length  submersed  In  liquid  nitrogen  upstream  to  the  discharge 

apparatus.  The  complete  apparatus  was  purged  for  at  least  12  hours  prior  to 

operation  with  a  slow  flow  of  99.993%  02  (<  10  p.p.m.  HgO) .  We  estimate 

residual  H20  density  from  inboard  leaks  to  be  less  than  1  p.p.m.  at  all  flow 

rates  used,  based  on  leak  rate  and  humidity,  but  quantitative  measurements 

should  be  made  in  the  Phase  II  experiments. 


Heterogeneous  Quenching  of  02(1z)  and  02(*a) 

«  1 

It  has  been  previously  establ1shedt■,  that  02(  z)  quenches  on  Pyrex  glass 

-2 

surfaces  with  a  quenching  efficiency,  y,  of  10  -  10  .  The  effective  rate 

constant  for  quenching  of  an  excited  species  flowing  In  a  tube  of  cylindrical 
cross  section  of  radius  r  is  given  by: 


rr2P 


2r,-l 
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Where  P,  D0  and  u  represent  the  overall  pressure,  the  diffusion  coefficient 
and  the  mean  molecular  speed,  respectively.  The  relative  contribution  of  each 
term  In  kHE  Is  summarized  In  Table  2  for  several  total  pressures  in  the  range 
0.1  -  10  Torr.  For  02(*a)  at  all  pressures,  the  second  term  in  Equation  2 
dominates  and  hence  the  actual  process  of  heterogeneous  quenching  at  the  wall 
Is  the  rate  limiting  step.  Using  the  values  in  Table  2,  kHE  <  0.01  sec"1  so 
that  the  heterogeneous  quenching  lifetime  is  greater  than  100  sec.  Since  the 
transit  time  of  a  given  volume  element  through  the  observation  vessel  ranged 
from  2  to  20  seconds,  the  depletion  of  09(^a)  by  heterogeneous  quenching  is 


clearly  negligible  In  the  present  experiments. 

For  O^1!),  heterogeneous  quenching  passes  through  two  regimes,  changing 
from  quenching  rate  limited  at  0.1  Torr  to  diffusion  rate  limited  at  10 
Torr.  The  present  experiments  were  performed  In  the  1-4  Torr  range,  where 
heterogeneous  quenching  Is  neither  strictly  diffusion  controlled  nor  quenching 
rate  controlled.  At  1.0  Torr  kHE  =  0.8  sec"1,  giving  an  effective  lifetime 
of  -  1.3  seconds  In  the  vessel.  At  1  Torr  the  pseudo-first  order  homogeneous 
rate  constant  for  quenching  of  02(1e)  by  ground  state  02  Is  1.3  sec-1.  At  4 
Torr  kHE  !  o.2  sec"1  while  the  rate  constant  for  quenching  by  ground  state 
Increases  to  5.2  sec-1.  With  a  transit  time  of  2  to  20  seconds  In  the  vessel, 
heterogeneous  quenching  of  02(1£)  Is  clearly  Important  at  pressures  of  -  1 
Torr  while  at  4  Torr  homogeneous  quenching  dominates.  In  the  limit  where  wall 
quenching  Is  diffusion  controlled,  the  density  profile  of  02(1e)  within  the 
tube  Is  expected  to  be  parabolic,  with  maximum  density  In  the  tube  center  and 
zero  density  at  the  wall.  We  have  measured  02(1c)  density  at  different 
distances  from  the  center  line  In  the  observation  vessel,  taking  Into  account 
the  diminishing  volume  viewed  as  the  detector  Is  moved  toward  the  edge  of  the 
tank.  The  results,  shown  In  Figures  2-4  show  that  at  P  *  0.4  and  0.8  Torr 
there  Is  no  clearcut  evidence  of  a  O^1!)  density  gradient  between  the  tube 
center  and  the  wall.  Nevertheless,  wall  quenching  of  02(1r)  Is  expected  to  be 
Important  as  pressures  approach  1  Torr.  The  results  could  be  explained 
If  y  «  2.5  x  10“3  In  Equation  (12). 


Flow  System  Characteristics 

The  time  base  used  for  the  monitoring  of  02(1a)  and  02(1e)  concentrations 
Is  determined  by  the  linear  flow  speed  of  02  in  the  observation  vessel  and  the 
distance  downstream  In  the  direction  of  flow.  This  simple  method  for  defining 


the  time  evaluation  of  the  composition  of  a  given  volume  element  of  gas  Is 
strictly  rigorous  only  In  the  case  of  “plug"  flow  where  the  diffusion  rate  In 
the  direction  of  flow  Is  much  smaller  than  the  linear  flow  rate. 

For  the  present  system,  we  have  attempted  to  Initially  create  plug-like 
flow  In  the  observation  chamber  by  Installation  of  polycarbonate  honeycomb 
flow  stralghteners  covered  on  one  side  with  80  micron  Teflon  screen  (see 
Figure  1).  To  observe  the  effect  of  these  devices  on  flow  characteristics  we 
rendered  the  flow  visible  by  the  addition  of  N£  to  the  O2  flow  upstream  of  the 
discharge.  The  ensuing  NO2  glow  Is  long  lived  and  easily  visible  with  the 
room  darkened.  We  were  able  to  determine  qualitatively  that  the  flow  stralght¬ 
eners  effectively  reduced  the  kinetic  energy  of  the  gas  flow  from  the  2"  lead- 
in  pipes,  converting  dynamic  pressure  to  static  pressure.  The  NO2  glow 
uniformly  filled  the  bell  shaped  reducer  leading  Into  the  observation  chamber. 
The  second  baffle  was  designed  to  further  reduce  the  flow  kinetic  energy. 

Without  the  baffles,  we  observed  turbulent  streaming  of  the  N02  glow  Into 
the  center  of  the  tank,  with  stagnation  occurring  at  the  tank  edges.  Although 
we  have  no  quantitative  measure  of  how  successfully  we  have  approached  plug 
flow  In  the  observation  vessel,  we  have  succeeded  in  at  least  uniformly 
filling  the  vessel  with  02  from  the  discharge,  as  substantiated  by  our  visible 
observations  and  by  the  Og^i)  concentration  profiles  shown  In  Figures  2-4. 
Measurements  with  a  hot  wire  anemometer  would  help  to  establish  the  degree  of 
flow  speed  uniformity. 

At  room  temperature  and  1  Torr,  the  diffusion  coefficient  for  02 
2  -1 

Is  -  150  cm  sec  .  Under  these  conditions  the  average  distance  a  molecule  of 
O2  will  diffuse  Is  -  25  cm.  At  linear  flow  speeds  of  5  -  15  cm/sec,  diffusion 
represents  a  significant  deviation  from  the  Idealized  method  of  extracting 
time  from  flow  speed  and  distance  measurements.  Since  diffusion  occurs  both 


In  the  forward  and  backware  direction,  the  net  effect  Is  to  reduce  concentra¬ 
tion  gradients  and  yield  apparent  species  decay  rates  which  are  smaller  than 
the  true  rates.  We  have  not  corrected  for  the  effects  of  diffusion  in  the 
present  experiment.  Adding  a  buffer  gas  diluent  should  reduce  the  diffusion 
coefficient  without  significantly  affecting  quenching  rates. 


C.  Kinetic  Analysis 

The  differential  rate  equation  governing  the  production  and  loss 
of  02(*a)  Is  given  as  follows: 

^Jr1  -- V3‘ll‘‘l  -2  V1»i2-kHE|lal  (13> 

♦  k^l'ni3!!  ♦ 

In  Equation  (13)  the  a)  loss  terms  represent  the  rate  of  quenching 
of  O^a)  by  ground  state  02,  wall  quenching,  and  energy  pooling,  while  the 
positive  terms  reflect  production  of  02(*a)  via  quenching  of  02(*i)  by 
02(3£)  and  by  02(1a). 

The  rate  equation  governing  production  and  loss  of  02(1e)  Is  given  by: 

=  -  k£l(3i)ll£)  -  kE2[1A)I1£]  -  kx[X)(l£| 

-k^l1!!  +  kpl^I2  (14) 

with  k£2  ■  k£2a  +  k£2b. 

1  1 
The  02(  £)  loss  terms  In  Equation  (13)  reflect  quenching  by  02(  £), 

quenching  by  O^a)  and  quenching  by  trace  Impurities,  l.e.  H20,  03,  0  etc. 

and  wall  quenching.  The  only  production  mechanism  for  02(1£)  Is  energy 

pooling.  Comparison  of  terms  In  Equations  (13)  and  (14)  reveal  the  assump- 


tlons  that  the  first  term  In  Equation  (11)  reflects  quenching  of  02(*z) 
to  02(*a),  while  quenching  by  impurities  leads  to  the  ground  state. 

In  the  present  system,  11a1/|3z1  <  0.08  and  [ 1/ I3zl  <  10“2,  so  that  the 
last  two  terms  in  Equation  (13)  may  be  neglected  unless  and  k£2a  are  very 
large.  Using  the  published  rate  constants  for  pooling  and  quenching  by  ground 
state  O2  (Table  1),  It  is  easily  shown  that  at  1  Torr  total  O2  pressure, 
quenching  of  O^1*)  by  ground  state  exceeds  loss  by  pooling  at  02(*a) 
fractions  lower  than  0.08.  Thus  at  02(*a)  conversion  fractions  significantly 
lower  than  this  value  (l.e.  <  0.005)  we  may  write  Equation  (13)  without 
serious  loss  of  accuracy  as: 


iUdtL  ■  -  ‘a  l3ml*i 
--ka  I1*' 


Equation  (15)  may  be  Integrated,  giving  a  pseudo-first  order  expression 
describing  the  decay  of  02(1a) : 


l1*!  =  l1*]^11 


As  shown  In  Section  B,  an  evaluation  of  the  magnitude  of  the  terms  in 
Equation  (14)  reveals  that  the  rate  of  heterogeneous  quenching  of  02(*z)  by 
the  vessel  wall  Is  of  the  same  order  of  magnitude  as  quenching  by  02(3z)  at 
pressures  below  1  Torr.  Similarly,  trace  quantities  of  H2O,  O3  and  0  are 
expected  to  be  significant  quenchers  of  02(^e)  at  densities  of  a  few  parts  per 
million.  Applying  the  steady  state  criterion  to  the  evolution  of  02(^z) 
density,  we  obtain  the  following  expression: 

-  -  k^zl^zl  -  kE2I1All1r]  -  kxlXlIl£l  -  kfc[h\  +  k  l1*!2  =  0  (17) 
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where  a  Is  an  arbitrary  constant.  Using  the  Initial  condition,  [ae1(0)  =  0, 
we  obtain: 


Til  (t)  = 
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Equation  (25)  predicts  that  02(*£)  density  will  build  up  to  a  maximum 

value  with  time  constant  determined  by  the  total  quenching  rate  of  O^r), 

+ 

then  decay  to  zero  with  time  constant  2kAl*  This  analysis  Is  similar  to 
those  given  by  O'Brien  and  Myers33  and  Heldner  et  al35  but  Includes  ral(t)  as 
a  quenching  species. 

The  effect  of  quenching  of  02(1i)  by  02(*a)  Is  expressed  In  the  time 

dependence  of  the  denominator  of  the  first  term  In  Equation  (25). 

If  kE2(A]Q  *  kg,  the  effective  rate  constant  for  decay  of  02(*i)  would 

*  *  1 

Initially  exceed  2  k^,  then  approach  2  k^  as  02(  a)  decayed  away. 

Although  not  attempted  In  the  present  Investigation,  a  trace  of  l£  may  be 
Introduced  Into  the  system  In  order  to  enhance  the  pooling  rate.  Kinetics 
studies34-3®  pertaining  to  operation  of  the  chemical  oxygen-iodine  laser  have 
led  to  the  elucidation  of  the  following  reaction  sequence: 
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The  rate  constant  for  Reaction  (28)  Is  reported^®  to  be  6  x  10^  times  greater 
than  that  for  Reaction  (1),  and  thus,  any  disadvantages  attending  the  addition 
of  Iodine  might  be  offset  by  the  enhanced  pooling  rate. 

In  the  COIL  laser,  energy  pooling  is  a  major  detriment.  Reaction  (1) 
depletes  O^a)  between  the  generator  and  the  laser,  while  Reaction  (28)  Is  a 
dominant  quenching  pathway  once  Iodine  Is  mixed  into  the  flow.  Under  these 
conditions,  a  steady  state  develops  described  by 


-  y1*]2  +  k28I1A]  11*1  -  (kQ  +  k£2I1A])[1r]  (29) 
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Since  the  equilibrium  between  I*  and  021*a]  is  established  very  rapidly: 


and  thus 
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When  I  Is  not  present  in  the  mixture.  Equation  (32)  of  course  reduces  to 
Equation  (18). 


0.  Results  and  Discussion 


We  have  fit  the  data  to  the  first  order  decay  of  02(*a)  and  02(*e) 
predicted  by  Equations  (16)  and  (25)  for  a  variety  of  flow  conditions.  A 
representative  decay  plot  for  02(*a)  and  02(*e)  Is  shown  In  Figure  5.  Within 
experimental  uncertainty  all  such  plots  exhibited  linearity,  with  the  ratio  of 

the  slopes  m(1z)/m(1A)  approximately  equal  to  2  as  predicted  by  Equation  (25). 

*  3 

The  dec*y  constants  (kAl  *  kAlrr])  derived  from  the  slopes  of  these  plots 
give  values  of  kAl  which  fall  Into  the  range  of  values  obtained  by  previous 
workers.  These  results  suggest  that  02(*a)  Is  being  quenched  largely  by 
02(^i)  (see  Table  1).  Another  Important  conclusion  obtained  from  these 
results  Is  that  we  are  unable  to  observe  the  density  maximum  predicted  by 
Equation  (25).  Not  even  at  the  highest  linear  flow  speeds  attempted  (100 
cm/sec)  were  we  able  to  resolve  the  rise  or  maximum  In  the  02(*e)  popula¬ 
tion.  This  could  be  because  the  maximum  occurred  In  the  upstream  portion  of 
the  apparatus  beyond  the  viewing  range  of  the  PMT.  A  direct  measurement  of 
the  02(*i)  risetime  would  yield  a  time  constant  t  which  reflected  the  total 
pseudo-first  order  rate  constant  for  02(*z)  destruction: 

■k-.t  . 

i  •  (Kq  +  k,i»0«  “  1  1  (33) 

In  order  to  assess  the  magnitude  the  rate  constant  for  quenching 
of  02(*i)  by  02(*a)  It  was  necessary  to  establish  the  magnitude  of  quenching 
by  other  species,  especially  those  whose  concentration  changes  with  time  (e.g. 
O3,  0).  In  the  absence  of  all  quenchers  except  ground  state  O2*  Equation  (19) 
predicts  that  the  steady  state  ratio  of  02(1t)  to  02(1a)  Is  given  by 


The  r.h.s.  of  Equation  (34)  thus  determines  the  upper  limit  of  the 
02(*e)  to  02(*a)  ratio.  A  quenching  species  whose  density  was  time 
Independent  would  suppress  I Js$  but  would  not  affect  the  rate  constant 
describing  Its  rate  of  decay  In  the  observation  vessel  at  long  times.  This 
rate  of  decay  Is  determined  by  k^.  A  quenching  species  whose  density 
decreased  with  time  would  similarly  suppress  |*e1,  but  the  magnitude  of  the 
effect  would  decrease  with  time. 

If  kl2  -  kEj,  the  resulting  plot  of  {*e!  vs.  time  would  be  expected  to 

exhibit  a  monotonic  decay  at  t  »  where  tm  Is  the  time  at  which  l1!) 

reaches  a  maximum.  In  this  regime  the  decay  constant  describing  02(*z) 

* 

disappearance  would  be  greater  than  2k  but  would  approach  this  value  as 

the  density  of  the  quencher  decreased  with  time. 

~  3  1 

Since  1 02 1  =  (  eJ  +  (  a  I  we  can  write  as  Equation  (18) 


(k2E  ‘  klE)  1  A|  +  k^10’1  +  Mxl  +  k' 
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which  can  be  rearranged  to 


=  k1  +  ks  l1*!"1 


klr  l°?l  +  k*  tXl  +  kHF 

where  k1  »  (k2j.  -  klE)/kp  and  k$  =  — - — ^ - — 
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A  plot  of  the  ratio  ( 1 o I /( I  vs.  is  expected  to  yield  a  straight  line 

with  slope  k$  and  Intercept  k^.  Thus  from  one  graph  we  may  extract  the 

1  3 

magnitude  of  quenching  of  02(  r)  by  the  wall  and  by  species  other  than  02(  e) 


from  the  slope  and  k^  from  the  Intercept.  Such  plots  obtained  under  several 
sets  of  conditions  are  shown  In  Figures  6-8.  As  these  figures  show  there  Is 
some  noise  In  the  plotted  data.  It  Is  suspected  that  the  major  source  of  this 
noise  Is  the  distortion  caused  by  the  nonunlform  glass  surface  of  the  Pyrex 
observation  vessel.  Visually,  the  vessel  surface  Is  slightly  rippled  and  Is 
far  from  good  optical  quality.  Furthermore  the  vessel  walls  become  thicker 
toward  the  ends,  and  some  attenuation  at  1.27  urn  was  noted  at  these 
positions.  For  this  reasons  data  were  taken  but  not  reported  at  positions 
where  attenuation  was  notice-able.  The  experimentally  determined  values  of  k^ 
are  listed  In  Table  3. 

Values  of  k^,  k$  and  their  uncertainty  limits  were  evaluated  using  a 
standard  linear  regression  statistics  package  (IMSL).  The  uncertainties 
represent  one  standard  deviation  at  the  (t  68X)  confidence  level.  Values 
of  k£2  calculated  from  using  Equation  (36)  are  giv^n  in  Table  3.  The 
magnitude  of  the  standard  errors  associated  with  the  Independent  measurements 
of  k£2  are  relatively  large,  however  averaging  all  independent  determinations, 
we  obtain  k^  *  6.5  ±  6.1  x  10“*®  cm'*  molec"*  $ec“*. 

Within  the  statistical  uncertainty  of  our  measurements,  we  conclude  that 
the  rate  constant  k^  for  quenching  of  by  Is  comparable  to 

krl*  the  rate  constant  f°r  Quenching  by  ground  state  oxygen.  Possible  sources 
of  error  In  the  estimation  of  k^  Include: 

1.  Uncertainty  In  02(*a)  and  Oj^e)  density  measurements  arising  from 
optical  variations  In  the  viewing  vessel. 

2.  Systematic  errors  due  to  quenching  of  02(*e)  by  residual  O3  or  oxygen 
atoms . 

3.  Errors  In  Interpretation  due  to  axial  diffusion  within  the 
measurement  vessel. 


Our  experimentally  determined  values  of  k$  obtained  via  Equation  (36) 
yield  values  of  k£^  which  agree  well  within  previously  determined  values 
of  k£j  ,  which  suggests  that  quenching  of  02(*z)  by  03  and  0  have  been 
effectively  eliminated  and  that  for  P  >  4  Torr,  heterogeneous  quenching  Is 
unimportant. 

The  evaluation  of  any  error  Introduced  by  axial  diffusion  of  components 
within  the  vessel  Is  difficult.  However,  since  the  measured  values  of  the 
rate  constants  for  quenching  of  (^(^z)  and  02(*a)  by  ground  state  oxygen 
(l.e.  kAl  and  k£^)  are  In  the  range  of  previous  determinations.  It  appears 
that  the  effect  Is  small.  The  fact  that  diffusion  distances  In  1  sec  at 
1  Torr  do  not  greatly  exceed  distances  traveled  by  the  flow,  would  tend  to 
support  this  conclusion,  however  further  Investigation  Is  needed  In  Phase  II. 

After  evaluating  these  sources  of  error,  we  have  concluded  that  the 
variable  viewing  efficiency  through  the  glass  vessel  Is  the  major  source  and 
that  the  resulting  errors  In  density  measurement,  though  systematic  In  fact, 
appear  random,  at  least  at  axial  positions  away  from  the  ends  of  the  vessel. 

In  the  absence  of  heterogeneous  quenching  and  quenching  by  trace 
Impurities, 


k£l[021 

k 


(37) 


A  measured  value  of  k$  greater  than  this  lower  limit  indicates  that  02(*£)  is 
probably  being  quenched  by  some  pathway  other  than  quenching  by  ground  state 
oxygen.  As  shown  In  Table  3,  the  derived  values  of  k£j  are  somewhat  higher 
than  the  literature  value  at  P  <  I  Torr,  but  at  P  >  4  Torr,  the  agreement  is 


very  good. 

We  believe  that  this  Is  the  first  singlet  oxygen  experiment  to  achieve 


02(  i)  densities  limited  only  by  homogeneous  quenching  by  ground  state  02. 
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The  number  densities  of  02(  £)  reported  here  (>  10  cm  )  appear  to  be  the 
highest  yet  achieved.  Other  Investigations  have  been  conducted  in  small 
diameter  flow  systems  where  wall  quenching  was  the  dominant  source  of  02(1£) 
deactivation. 

E.  Suimnary 

We  have  constructed  and  successfully  demonstrated  a  discharge-flow 
apparatus  for  02(1a)  and  02(*i)  kinetics  studies.  Some  important  results 
obtained  Include: 

1)  Verification  of  the  rate  constant  for  quenching  of  02(1a)  by  ground 
state  oxygen. 

2)  Attainment  of  02(1z)  concentrations  limited  predominantly  by  the  rate 
of  quenching  of  02(1r)  by  ground  state  02.  We  believe  that  this  Is  the 
first  discharge-flow  experiment  to  produce  02(1r)  under  conditions  where 
wall  quenching  or  quenching  by  Impurities  were  not  the  dominant  factors 
limiting  the  steady  state  density  of  02(1s). 

3)  Estimation  of  k^,  the  rate  constant  for  quenching  of  02(1z)  by 
02(*a).  Although  the  value  obtained  for  kz2  is  somewhat  uncertain  due  to 
noisy  data.  Its  magnitude  (6.5  ±  6.1  x  10-16  cm3  molec-*1  sec-1)  is  not 
extremely  large  compared  to  k£^,  the  rate  constant  for  quenching  by  ground 
state  02.  We  find  these  preliminary  results  very  encouraging  as  they 
pertain  to  the  feasibility  of  the  Red  SOX  laser. 
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Table  3.  Rate  constants  derived  froa  experimental  data 
(see  Figures  6-8).  Values  for  k  .  and  k  « 
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Figure  2.  Relative  0,(  E)  density  measured  at  various  distances 
from  tube  axis.  Plottsd  on  the  horizontal  axis  Is 
distance  in  the  direction  of  the  flow  axis. 
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Figure  5.  Tine  dependence  of  emission  signal  fron  0.(1A) 
(opes  dieaonds)  end  02 (i) (filled  diamonds) . 


